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Magnetism  in  the  orthorhombic  metal  CaFe4As3  was  examined  through  neutron  diffraction  for 
powder  and  single  crystalline  samples.  Incommensurate  (qm  ~  (0.37  —  0.39)  xb*)  and  predominantly 
longitudinally  (||  6)  modulated  order  develops  through  a  2nd  order  phase  transition  at  Tn  =  89.63(6) 

K  with  a  3D  Heisenberg-like  critical  exponent  /?  =  0.365(6).  A  1st  order  transition  at  T2  =  25.6(9) 

K  is  associated  with  the  development  of  a  transverse  component,  locking  qm  to  0.375(2)6’ ,  and 
increasing  the  moments  from  2.1(1)  to  2.2(3)  ps  for  Fe2+  and  from  1.3(3)  to  2.4(4)  pb  for  Fe+.  The 
ab-initio  Fermi  surface  is  consistent  with  a  nesting  instability  in  cross-linked  FeAs  strips. 

PACS  numbers:  75.30. Fv,  64.70.Rh,  61.05.F-,  71.18.+y 


The  recent  discovery  of  superconductivity  (SC)  in  iron 
pnictides  [lj  has  opened  a  new  chapter  in  SC  research. 
In  contrast  to  the  cuprates,  which  display  SC  near  a  cor¬ 
related  insulating  state,  iron  SC  appears  when  doping  or 
applying  pressure  to  metallic  parent  compounds  with  a 
spin  density  wave  (SDW).  A  square,  puckered  iron  chal- 
gogenide  or  iron  pnictide  plane  is  the  common  theme 
in  this  family  of  materials,  which  so  far  has  been  real¬ 
ized  in  the  ZrCuSiAs  (1111),  ThCr2Si2  (122),  anti-PbO 
(11)  and  Cu2Sb  (111)  structures.  Chemically  related 
CaFe4As3  ii  contains  interpenetrating  FeAs  strips, 
which  -as  spin-ladders  in  copper  oxides-  may  provide 
unique  insight  into  electronic  correlations  of  the  square 
lattice  plane. 

Similar  to  the  parent  compounds  of  the  1111-  and  122- 
type  superconductors,  CaFe4As3  is  not  superconducting 
down  to  1.8  K  [3|.  The  novel  framework  structure  is 
composed  of  shared  FeAs  tetrahedra  with  Ca  atoms  in 
the  channels  so  defined  (Fig.  1(a)).  The  structure  is 
reminiscent  of  porous  manganese  oxides  0].  Extending 
along  the  6-axis  are  FeAs  strips  where  4-fold  As  coordi¬ 
nated  Fe  (Fel-Fe3  in  Table  I)  occupy  an  approximately 
square  lattice  of  finite  width.  Fe  atoms  that  link  strips 
(Fe4  in  Table  I)  on  the  other  hand  form  FeAss  pyramids. 
Mossbauer  spectra  [3]  indicate  different  electronic  states 
on  4-  and  5-fold  coordinated  sites  consistent  with  formal 
valences  of  Fe2+  and  Fe1+,  respectively. 

Two  phase  transitions  were  detected  through  thermal 
anomalies  in  bulk  properties  ii-  While  only  the  up¬ 
per  transition  at  Tn  =  89.63(6)  K  is  marked  by  a  sig¬ 
nificant  specific  heat  anomaly,  the  lower  transition  at 
T2  =  25.6(9)  K  yields  a  strong  reduction  in  resistivity 
(p).  The  magnetization  is  anisotropic  with  larger  val¬ 
ues  along  strips  (||  b).  The  large  Sommerfeld  constant 


(7  =  0.02  J/mol-Fe/K2)  indicates  relatively  strong  elec¬ 
tronic  correlations. 

To  elucidate  the  magnetic  properties  of  CaFe4As3,  we 


FIG.  1:  (color  online).  Crystal  (a,b)  and  magnetic  structures 
of  CaFe4As3  at  (c,d)  T  =  50  K  and  (e,f)  5  K  in  ac  plane  (left) 
and  along  6-axis  (right).  Solid  squares  indicate  a  unit  cell. 
The  vector  d  (a  ■  d  =  cos  30°)  is  used  for  right  panels. 
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TABLE  I:  Atomic  positions  within  Pnma  of  CaFe4As3  at  RT 
determined  by  Rietveld  analysis  (y2  =  2.2).  Lattice  constants 
are  a  =  11.9062(2)  A,  b  =  3.74477(6)  A,  and  c  =  11.6131(2) 
A.  Isotropic  Debye-Waller  factors  (t/iso)  are  employed. 


Atom 

Site 

X 

y 

z 

Uiso  (A2) 

Asl 

4  c 

0.1364(4) 

1/4 

0.0778(4) 

0.0149(12) 

As2 

4c 

0.3839(4) 

1/4 

0.2636(4) 

0.0093(13) 

As3 

4c 

0.0998(4) 

3/4 

0.3797(3) 

0.0050(13) 

Fel 

4c 

0.1808(2) 

1/4 

0.2763(2) 

0.0055(9) 

Fe2 

4c 

0.1941(3) 

1/4 

0.8759(3) 

0.0095(8) 

Fe3 

4c 

0.4343(3) 

1/4 

0.4635(3) 

0.0064(9) 

Fe4 

4c 

0.4770(3) 

1/4 

0.6864(3) 

0.0103(10) 

Ca 

4c 

0.3376(4) 

3/4 

0.0686(6) 

0.0115(16) 

carried  out  neutron  scattering  measurements  on  powder 
and  single  crystalline  samples.  We  show  that  the  up¬ 
per  transition  is  to  a  longitudinal  incommensurate  (IC) 
SDW,  while  the  lower  transition  is  associated  with  the  de¬ 
velopment  of  a  transverse  component  and  lock-in  to  com¬ 
mensurate  order.  A  powder  sample  and  a  single  crystal 
of  CaFe4As3  were  prepared  according  to  Ref.  Neu¬ 
tron  powder  diffraction  data  were  collected  on  BT1  and 
BT9  at  NIST.  Single  crystal  experiments  were  performed 
on  BT9  and  on  the  newly  developed  Multi  Axis  Crystal 
Spectrometer.  The  error  bars  indicate  ±  the  standard 
deviation.  We  used  group  theoretical  analysis  to  identify 
magnetic  structures  that  are  allowed  by  symmetry  and 
consistent  with  the  diffraction  data. 

We  performed  neutron  powder  diffraction  measure¬ 
ments  at  T  =  5  K,  50  K,  100  K  and  room  temperature 
(RT).  As  in  previous  reports  Hi  ,  the  crystal  structure 
is  accounted  for  by  space  group  Pnma.  Table  U  sum¬ 
marizes  the  RT  structural  parameters.  Additional  peaks 
associated  with  magnetic  ordering  were  observed  at  50  K 
and  5  K.  Figure  2  shows  a  portion  of  the  Rietveld  refined 
diffraction  pattern,  which  is  consistent  with  the  single 
crystal  data.  All  the  magnetic  peak  positions  can  be  ac¬ 
counted  for  by  a  propagation  wave  vector,  qm  =  (OkO). 
The  strongest  magnetic  peak  actually  consists  of  over¬ 
lapping  (1 A0)  and  (Ofcl)  peaks,  where  k  =  0.386(2)  and 
0.375(2)  at  50  K  and  5  K,  respectively,  indicating  qm  is 
temperature  dependent. 

To  determine  the  temperature  dependence  of  qm,  we 
used  single  crystal  diffraction  in  the  (0 kl)  plane  on  a 
0.5  x  8  x  0.5  mm3  sample.  Figure  3  shows  the  In¬ 
dependence  of  the  (Ofcl)  magnetic  peak.  The  peak  ap¬ 
pears  for  T  w  90  K  along  with  an  anomaly  in  p(T)  |2|. 
The  temperature  dependent  integrated  intensity  (/)  can 
be  described  as  I  oc  f2/3,  where  t  =  (Tn  —  T)/Tn  is  the 
reduced  temperature  and  Tn  =  89.63(6)  K  (right  inset  to 
Fig.  3(b)).  The  inferred  critical  exponent  (3  =  0.365(6) 
is  consistent  with  the  3D  Heisenberg  model,  0.36  [§]. 

Below  Tn,  the  peak-center  k  appears  to  increase  con¬ 
tinuously  from  0.37  to  0.39  on  cooling  (Fig.  3(a)),  which 
indicates  an  incommensurately  modulated  structure.  / 


and  k  then  change  abruptly  near  25  K.  We  associate  this 
anomaly  with  a  first  order  phase  transition  for  the  follow¬ 
ing  reasons:  (i)  There  is  considerable  thermal  hysteresis 
(AT  ~  2  K)  of  I  and  p  (Fig.  3(c))  and  (ii)  diffraction 
peaks  associated  with  the  two  phases  coexist  close  to  T2 
(left  insets  to  Fig.  3(b)).  Since  no  structural  phase  tran¬ 
sition  was  detected  at  T2,  this  appears  to  be  an  intrinsic 
first  order  magnetic  phase  transition.  This  is  different 
from  observations  in  122  systems,  where  a  first  order 
structural  phase  transition  precedes  @  or  coincides  [lj 
with  a  magnetic  transition.  While  the  structural  phase 
transition  in  122  systems  breaks  tetragonal  symmetry 
and  relieves  frustrating  magnetic  interactions  such  a  sce¬ 
nario  is  not  expected  in  orthorhombic  CaFe4As3. 

Below  T2,  k  is  T-independent  and  fixed  at  3/8  within 
experimental  accuracy.  This  indicates  an  IC  to  commen¬ 
surate  phase  transition.  Such  transitions  are  typically 
found  to  be  first  order  because  they  involve  the  loss  of  a 
soft  phase  degree  of  freedom  that  exists  in  the  IC  state 
but  is  lost  as  the  spin  structure  becomes  commensurate 
and  pinned  to  the  lattice  in  the  low-T  state.  The  tran¬ 
sition  also  involves  an  abrupt  drop  in  p  along  the  fo-axis 
which  is  consistent  with  loss  of  scattering  mechanisms 
associated  with  the  phase  degree  of  freedom  in  the  IC 
state.  The  transition  is  only  visible  as  a  weak  inflexion 
in  specific  heat  data  ,2|,  suggesting  that  the  dominant 
modification  at  T2  is  a  change  in  magnetic  periodicity, 
which  permits  phase-locking  the  SDW  to  the  lattice. 

For  insights  into  the  origin  of  incommensurability,  we 
calculated  the  electronic  band  structure  using  the  den¬ 
sity  functional  in  an  all  electron  full  potential  linearized 
augmented  planewave  method  with  local  orbitals.  We 
used  the  WIEN2K  [7j  code  with  a  generalized  gradient 
approximation  and  spin  orbit  interactions.  The  corre¬ 
sponding  Fermi  surface  is  shown  in  Fig.  4  (a).  While  the 
local  chemistry  of  Fe  is  similar  to  that  of  other  pnictides, 


QW1) 


FIG.  2:  (color  online).  A  portion  of  the  CaFe4As3  neutron 
diffraction  pattern  at  100  K,  50  K  and  5  K  with  Rietveld 
refinement  (solid  lines).  The  black  and  green  bars  with  repre¬ 
sentative  (hkl)  indices  show  the  calculated  positions  of  nuclear 
and  magnetic  peaks,  respectively. 
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FIG.  3:  (color  online),  (a)  Color  plot  of  the  neutron  scat¬ 
tering  intensity  versus  temperature  and  wave  vector  transfer 
along  the  (Ofcl)  direction.  The  temperature  dependence  of 
the  inferred  integrated  intensity  and  peak  position  are  shown 
in  frames  (b)  and  (c)  respectively.  The  vertical  bar  in  (c)  rep¬ 
resents  a  systematic  instrumental  uncertainty.  The  insets  to 

(b)  show  warming  intensity  against  the  reduced  temperature 
near  Tn  and  coexistence  of  two  peaks  near  T 2.  The  insets  to 

(c)  emphasize  thermal  hysteresis  of  I  and  p. 


the  fermiology  is  distinct.  Consistent  with  the  strip-like 
structure  and  the  quasi-ID  resistivity,  the  dominant  fea¬ 
tures  are  sheets  parallel  to  the  6-axis.  As  in  the  122  pnic- 
tides,  there  are  electron  and  hole  tubes  (Fig.  4(b))  near 
the  corners  of  the  Brillouin  zone  and  T  point  respectively. 
The  wave  vector  connecting  the  tubes,  however,  has  no 
relation  to  magnetic  ordering.  The  dominant  nesting  in¬ 
stability  appears  to  be  associated  with  the  planar  Fermi 
sheets  which  are  displaced  from  each  other  by  distances 
close  to  3b*/8.  This  indicates  an  itinerant  description  of 
magnetism  and  Fermi  surface  nesting  may  be  relevant. 

We  now  turn  to  the  detailed  nature  of  the  two  magnetic 
phases  at  50  K  and  5  K  using  a  single  crystal  oriented 
in  (hkO)  and  (0 kl)  planes.  The  observed  intensities  were 
compared  to  the  intensities  calculated  for  various  irre¬ 
ducible  representations  (irreps)  (s|  in  Figs.  5.  The  50  K 


FIG.  4:  (color  online),  (a)  Fermi  surface  of  CaFe4As3.  Unlike 
in  the  122  family,  plane-like  Fermi  surfaces  are  dominant,  (b) 
Viewing  along  kc  reveals  tube-like  pockets  as  well  as  plane-like 
Fermi  surfaces.  The  nesting  vector  k  =  3/8  is  indicated. 


and  5  K  spin  structures  corresponding  to  the  best  fits  are 
shown  in  Figs.  l(c,d)  and  l(e,f),  respectively.  At  50  K, 
the  data  are  consistent  with  magnetic  moments  (||  6)  in 
irrep  74.  Components  within  the  ac  plane  are  allowed  by 
symmetry  but  not  supported  by  the  data  (0.04(13)  /xb). 
The  wave  vector  at  this  temperature  is  (0,  0.386(2),  0) 
and  the  spin  structure  may  be  described  as  a  longitudi¬ 
nal  SDW.  The  best  fit  moment  sizes  are  2.1(1)  piB  and 
1.3(3)  hb  per  Fe2+  and  Fe+  sites,  respectively. 


At  5  K,  below  the  first  order  transition,  qm  is  indistin¬ 
guishable  from  (0,  3/8, 0),  and  the  refinement  indicates  an 
ac  plane  spin  component  described  by  irrep  t\  appears 
below  T2.  A  full  four-circle  experiment  to  definitively  de¬ 
termine  the  ordered  spin  structure  awaits  availability  of 
larger  single  crystals.  The  sequence  of  two  phase  tran¬ 
sitions  first  to  longitudinal  IC  order  followed  by  the  ap¬ 
pearance  of  a  transverse  component  with  a  different  irrep 
also  occurs  in  TbMn03  [f|  and  Ni3V20s  40],  where  it 
can  be  accounted  for  by  a  combination  of  competing  ex¬ 
change  interactions  and  easy  axis  anisotropy  as  described 
by  Ref.  11|.  A  longitudinal  to  transverse  SDW  transi¬ 
tion  is  also  often  found  in  itinerant  systems,  Cr  being  a 
familiar  example  0.  As  indicated  by  the  enhanced  low- 
T  intensity  (Fig.  3(b)),  the  moment  size  increases  below 
T2  and  the  analysis  yields  2.4(4)  /xb  per  Fe+  and  2.2(3) 
/xb  per  Fe2+ .  By  far  the  largest  increase  occurs  for  the  5- 
fold  coordinated  Fe+  site  that  links  FeAs  strips.  Overall, 
the  ordered  moment  is  substantially  larger  than  0.4-0. 9 
/xs/Fe  observed  in  the  1111  and  0. 8-1.0  /xs/Fe  in  the  122 
systems  [lj,  suggesting  stronger  electron  correlations. 


For  a  theory  of  magnetism  in  CaFe4  As3 ,  we  performed 
first  principles  total  energy  calculations.  The  antiferro¬ 
magnetic  (AF)  state  identified  experimentally  was  found 
to  have  a  lower  energy  than  the  ferromagnetic  and  para¬ 
magnetic  states.  Analysis  of  the  corresponding  staggered 
magnetization  reveals  that  while  the  fundamental  peri¬ 
odicity  prescribed  by  the  supercell  is  1/8,  the  dominant 
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FIG.  5:  (color  online).  Observed  ( x  axis)  and  calculated  ( y 
axis)  values  for  the  squared  magnetic  structure  factors  within 
representative  irreps  at  (a)  T  =  50  K  and  (b)  5  K. 


harmonic  is  3/8  in  good  agreement  with  the  experiments. 
While  the  antiferromagnetism  results  in  a  calculated  re¬ 
duction  of  the  density  of  states  at  the  Fermi  level  by  a 
factor  of  ~  2,  LDA  predicts  the  low-T  magnetic  state 
remains  metallic  as  observed  experimentally. 

The  structural  similarity  between  the  FeAs  strips  in 
CaFe4As3  and  FeAs  layers  in  the  1111,  122,  11  systems 
suggests  analogies  between  their  magnetic  structures.  We 
shall  describe  these  in  the  2D  centered  square  unit  cell 
that  holds  two  Fe  atoms  and  has  a  lattice  parameter  a/2 
times  the  Fe-Fe  spacing.  In  an  FeAs  strip  description  of 
CaFe4As3,  the  6-axis  corresponds  to  a  lattice  parameter 
in  that  cell.  The  most  common  spin  structure,  AF2,  has 
a  wave  vector  (0.5, 0.5)  and  corresponds  to  M-point  Fermi 
surface  nesting.  AF2  features  nearest  neighboring  (nn) 
spins  that  are  AF  correlated  in  one  direction  and  ferro- 
magnetically  correlated  in  the  perpendicular  direction.  It 
is  called  striped  with  stripes  extending  along  nn  bonds 
as  favored  by  next  nn  exchange  interactions.  AF2  is  ob¬ 
served  in  the  parent  compounds  of  the  122  and  1111  ma¬ 
terials  and  is  associated  with  the  spin  fluctuations  from 
which  a  spin  resonance  is  formed  in  the  superconducting 
state.  AF2  is  to  be  distinguished  from  the  conventional 
nn  AF,  which  we  denote  AF1  with  a  wave  vector  (1,0)  in 
the  centered  cell.  AF1  is  also  a  striped  structure  but  the 
stripes  extend  along  next  nn  bond  directions.  Favored  by 
nn  exchange  interactions,  AF1  is  consistently  observed  in 
the  parent  compounds  to  cuprate  superconductors  but 
not  so  far  in  any  compound  related  to  iron  pnictides.  Fi¬ 
nally,  AF3  has  a  wave  vector  (0.5,0)  and  is  typically  called 
double  striped,  stripes  extending  again  along  nn  bonds 
in  a  structure  that  is  favored  by  third  nn  exchange  inter¬ 
actions.  AF3  and  IC  variants  of  it  (6,0)  with  5  =  0.3— 0.5 
are  observed  in  Fei+yTei_j,Sex  for  small  x  and  finite  y , 
respectively  JT^].  AF3  is  analogous  to  the  structure  ob¬ 
served  in  CaFe4As3  where  5  =  3/8  at  low-T.  Indeed 
examining  Fig.  1(d)  one  can  appreciate  the  near  double 
striped  nature  of  the  structure.  A  similar  periodicity  is 
also  reported  for  structurally  related  Fei.uiTe  fl3|. 


It  has  been  pointed  out  that  thicker  FeAs  planes  are 
associated  with  the  AF3  structure  while  the  AF2  struc¬ 
ture  occurs  for  thinner  layers.  The  layer  thickness  is  3.6 
A  in  Fei.o76Te  [l3}  with  AF3  structure  and  2. 7-2. 8  A  in 
the  122  systems  with  the  AF2  structure.  LDA  calcula¬ 
tions  have  suggested  a  critical  layer  thickness  of  3.4  A  for 
the  transition  from  AF2  to  AF3  related  to  the  effect  of 
changing  bond  angles  on  second  and  third  nn  exchange 
interactions  0-  In  this  context,  it  is  interesting  to  note 
that  the  average  layer  thickness  in  CaFe4As3  ranges  from 
3. 0-3. 3  A;  substantially  exceeding  that  for  the  122  and 
1111  structures,  though  not  quite  reaching  the  inferred 
threshold  for  the  AF3  structure.  To  assess  the  relevance 
of  competing  exchange  and  local  moment  magnetism  ver¬ 
sus  Fermi  surface  nesting  and  an  itinerant  description, 
doping  studies  should  be  useful  as  in  metallic  Cr  [12].  A 
shift  in  the  wave  vector  with  carrier  doping  would  indi¬ 
cate  itinerant  magnetism  controlled  by  nesting. 


In  summary  we  have  examined  FeAs  magnetism  in  a 
unique  strip  geometry  realized  in  CaFe4As3-  A  longitu¬ 
dinal  SDW  with  IC  modulation  similar  to  observations 
in  FeAs  and  Fei+yTe  occurs  below  Tn  =  89.63(6)  K.  A 
first  order  transition  at  T2  =  25.6(9)  K  yields  a  trans¬ 
verse  component  and  locks  the  SDW  to  the  lattice  with 
a  commensurate  wave  vector  (3b* /8).  The  calculated 
Fermi  surface  features  sheets  separated  by  a  wave  vector 
similar  to  that  of  the  SDW,  suggesting  that  Fermi  sur¬ 
face  nesting  plays  a  role.  Alternatively  the  IC  state  could 
result  from  competing  second  and  third  nn  exchange  in¬ 
teractions  in  a  localized  spin  picture.  Dramatic  changes 
in  p  through  the  lower-T  phase  transition  indicate  prox¬ 
imity  to  3 d  electron  localization,  a  condition  that  appears 
favorable  to  high-Tc  SC. 
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